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Abstract 

The generalized defect D(J7) of a unitary N x N matrix U with no zero entries 
is the dimension of the real space of directions, moving into which from U we do 
not disturb the moduli \Ui j\ as well as the Gram matrix U*U in the first order. 
Then the defect d(U) is equal to D(£7) — (2N — 1), that is the generalized defect 
diminished by the dimension of the manifold {D r UD c : D r , D c unitary diagonal}. 
Calculation of d(U) involves calculating the dimension of the space in M, N spanned 
by a certain set of vectors associated with U. We split this space into a direct sum, 
assuming that U is a Kronecker product of unitary matrices, thus making it easier 
to perform calculations numerically. Basing on this, we give a lower bound on 
D(J7) (equivalently d(£7)), supposing it is achieved for most unitaries with a fixed 
Kronecker product structure. Also supermultiplicativity of D(C/) with respect to 
Kronecker subproducts of U is shown. 

Keywords: unitary matrix, kronecker product, tensor product, complex Hadamard 
matrix, doubly stochastic matrix, critical point 

MSC-class: 51F25, 47A80, 15B34, 15B51, 58A05 

1 Introduction 

The paper is a sequel to [T] where the defect of a unitary matrix was introduced and 
investigated in a general case. 

Let U G U, where by U we shall denote the set of all unitaries N x N . Consider the 
problem: find all directions (complex matrices) guaranteeing no first order change of \Uij\ 
for all i,j when moving from U in such a direction, and at the same time no first order 
disturbance of unitarity of U. The first requirement implies that these directions are of 
the form \RoJJ where o is the Hadamard product and R is a real matrix. The second 
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means that the directions in question must belong to the space tangent, at U, to the N 2 
dimensional manifold U. This tangent space is the set {EU : E antihermitian} or, if one 
likes, {UF : F antihermitian}. 

Note that the set of directions in question forms a real linear space and it is overpa- 
rameterized by R in the case of U having zeros among its entries. If U^j = then the 

rp 

equality i(R + aei(ej) )oU = iRoU, with being the standard basis column vectors, 
exposes the fact that two different i?'s can be used to parameterize any single direction. 
Having this in mind we thus concentrate on the space of R's parameterizing the space of 
directions in question: 

{R : \RoJJ = EU for some antihermitian E} (1) 

and ask about its dimension. We define the defect of U as this dimension diminished by 
27V- 1. 

If U has no zero entries, the defect plus 2iV — 1 is precisely the dimension of the space of 
directions in question. The reduction by 2iV — 1 is motivated by applications of the notion 
of defect. For most unitaries, especially those with no zero entries, the 2iV directions: 

\ei(efoU = idiag(ei)-£/ for i = 1..N, (2) 
i e(e_j) o U = U ■ i diag(ej) for j = 1..N, 
where e = ([1 1 1 . . . 1]) T 

which are also directions in question, span the 2iV — 1 dimensional real space tangent to 
the phasing manifold {D r UD c : D r , D c unitary diagonal} generated by left and right 
multiplication of U by unitary diagonals. The manifold consists of what we regard as 
unitaries equivalent to U, with the property of having moduli of entries identical with 
those of U. Usually the directions (J2]), leading from U to these equivalent unitaries, are 
uninteresting. 

Other directions in question may lead to unitaries inequivalent to U, but with the same 
pattern of moduli, provided they exist in a neighbourhood of U. A special case of this 
property is the existence of a smooth family of inequivalent complex Hadamard matrices 
stemming from a given complex Hadamard matrix H. Here by a complex Hadamard 
matrix we understand, unlike in combinatorics community, a unitary with all moduli 
equal to 1/ V~N, although matrices rescaled to have unimodular entries are more often 
considered in this context. As far as such complex Hadamard matrices are concerned, 
they have many applications especially in quantum information theory. We recommend 
the reader to consult the online catalogue of complex Hadamard matrices available at [3J. 

The smallest matrix size for which the known collection of inequivalent complex Hadamards 
is presumably incomplete is 6. Current investigations are concentrated on the search for 
a 4-dimensional smooth family of inequivalent complex Hadamards. This search is partly 
motivated by the fact that the value of the defect calculated for most known 6x6 complex 
Hadamards, including the unitary Fourier matrix Fq, is equal to 4. This is because the 
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defect of a complex Hadamard H provides an upper bound for the dimension of a smooth 
family of inequivalent complex Hadamards stemming from H. We suggest the reader to 
throw a look into our previous work [1] to find more details about this application of the 
defect. There he can also find that the defect of unitary U (complex Hadamard H) being 
equal to zero implies that there are no unitaries (complex Hadamards) inequivalent to 
U (H) in its neighbourhood. At this point it is worth to recall that analogous (see pQ) 
result was obtained by Remus Nicoara in his work [I], under the name of 'span condition', 
but in a different context of commuting squares of algebras, in which complex Hadamard 
matrices are used to define commuting squares of matrix algebras. 

This above mentioned value 4, as we now understand, causes one to suppose that all these 
6x6 complex Hadamards (or complex Hadamards equivalent to these) with such a defect 
are connected by a single 4-dimensional family. The old 2-dimensional family stemming 
from the Fourier F$ was reported, (among other even earlier sources [5], [6]) already in our 
catalogue [7J [3] . Other 1-dimensional families were constructed in [8j [9] , 2-dimensional 
in pilj IUJ, 3 dimensional in [TJ], and the first attempt to build a 4-dimensional family of 
6x6 inequivalent Hadamards has to be attributed to Ferenc Szollosi due to his succesful 
work [13], inspired also by some numerical evidence of the existence of such a family 
contained in \14 



From the above collection of results we conclude that it is of value to study the notion 
of defect, to search for effective ways to calculate it or estimate it for unitary matrices of 
a large size. 

In the following we employ another, more expedient characterization of the defect. To 
this end let us start with rewriting the definition flTJ), 



\RoJJ = EU for some antihermitian E 
(iR o U)(U)* is antihermitian 

u(Rouy = (Rou)(uy , 



(3) 



where * denotes the Hermitian conjugate. 
This further leads to 

V 1 < i < j < N [U],, [diag([%)] {[U] j>: ) = [f/]. : [diagQijyj [[U]^ 

t (4) 

V 1 < i < j < N [f/] v [diag([%)] ([£%)* = \U] t . [diag([i?] v )] ([[/] ^ 

where by [X] k . the k-th row of X is denoted, we adopt MATLAB notation for rows and 
columns. 

The above can be rewritten as 

Vl<Kj<iV tri^M; 



(id) 
U 



(5) 
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where 



... o ••• 



[U\i, ° W; 



M, 



[hi) 



A 



U 



i • S, 



(6) 



" PI, ° [U] 3 , 





has only its z-th and j-th rows potentially filled with nonzeros. Here A^' 3) and S } 
denote the real and imaginary part of M$ J) . 

Let us convert (jSJ) into purely real conditions: 

T 



V 1 < i < j < N tr R^A^ = and tr r(S ( 



(id) 

u 



(7) 



The form trA(B) is a Hilbert-Schmidt type inner product on the space of real N x N 
matrices, hence (J7|) amounts to orthogonality conditions in this space. Consequently, and 
in accordance with our initial definition the defect of U can be calculated as 

dim ({R : iRoU = EU for some antihermitian E}) — (2N — 1) 



dim ^span ^ 



.4 



(M) q(i,j) 

U ' °u 



: 1 < % < j < N 



(2N - 1) (8) 



where () ± is used to indicate the orthogonal complement of the space M.jj spanned by 



matrices Ay , Sy^, i.e. 



where, 



M 



v 



def 



span (Mi 



(9) 



M 



u 



def 



A% j) : 1 < % < j < n\ U {S { v' j) : 1 < i < j < N 



(10) 



The properties of orthogonal complements allow now for the following alternative defini- 
tion of the defect: 



Definition 1.1 The defect of a unitary N x N matrix U is the number 

d(U) = dim ((Mi/)- 1 ) - (2JV-1) 
= (N-l) 2 - dim (M^) . 



11^ 
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It was shown by Karabegov [2] that the set of those U for which the defect is greater 
than is of the Haar measure zero within the set of all unitaries N x N. In this article 
we give a lower bound on the defect of a Kronecker product of unitary matrices: 

U = U {1) <g> U {2) ® . . . ® U (r) , (12) 

where the unitaries are of size ri\ x n%, x 712, n r x n r , respectively. The bound 
computed according to our formulas agrees with the defect obtained numerically for 
Kronecker products of random unitary matrices drawn according to the Haar measure 
on the unitary group. 

However, we do not have here a result concerning measures, analogous to that of Karabegov 
Namely, that for unitary matrices lying beyond a certain set of the Haar measure 
zero, within the set of all unitary x matrices, the defect of the resulting Kronecker 
product ffT2]) is equal to our bound. 

Note that the above Kronecker product structure forms a key tool in studying quantum 
composite systems. For instance, any local unitary operation acting on a system com- 
posed of r particles has the form ( lT2"j) . Thus our results on the defect describe algebraic 
properties of generic local unitary operations, which are also called quantum gates. 



2 Preparation 



Let us introduce vector indices into Kronecker product (1121) . 



where 



Ph,M,^M = U ^ » ( 13 ) 



i = (ii - 1) JJ rifc + (« 2 — 1) JX n ft + ••• (v-i-l)"r + ir (14) 

k=2 k=3 

r r 

3 = Ul- l )\\_ n k + C72-l)IJw* + ••• (jV-l-l)^r + jr (15) 
k=2 k=3 

and 

[%..v],b,,V] = • P i2 \ h • ••• • [ Ui % r , jr (16) 

In the following, if not stated otherwise, i, j will correspond to the vector indices [ii..i r ], 
b'i-Jr]- 

Let us try for a while to understand vector indices better. The relation between ordinary 
and vector indices could be defined in a different, but completely equivalent way. For 
two distinct elements in {l..ni} x . . . x {l..n r } let us introduce a relation: 

[z'i..v] < bWr] <=^> h=ji, ■ ■ ■ , i P = j P , ip+i < jp+i for some pe{l,..,r-l} . 

(17) 
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Then the relation between ordinary indices and vector indices we choose to be the only 
bijection <fi from {1, . . . , YYi=i n i\ ' m ^° X ... X {l..n r } (for which we write 

4>{i) = [ii..i r ]), which satisfies: i < j ==>- <f)(i) < 

Later in this article we talk about subrows and submatrices, built from elementary objects 
(entries of a row for a subrow, rows or columns of a matrix for a submatrix) indexed by 
those i, whose k-th subindex i k is fixed at some value: i E A = {i : ik = a}. Any i G A is 
therefore determined by an element of x . . . x {1.. n k ^i} x {1. .rik+i} x . . . x {l..n r }, 

namely the reduced vector index [ii..ik-i, i k +i..i r ] being a subvector of [i\..i k -\, a, i k+ \..i r ] 
corresponding to i. We can introduce the relation < for the reduced vector indices, and 
then we note that: 

[ii..ik-i,ik+i~ir] < \ji~jk-ijk+i~jr] ^=^" [ii..i k _ 1 ,a,i k+1 ..i r ] < \ji..jk-i,oc,j k +i-jr] ■ 

(18) 

Ordering the reduced vector indices induces, again, their relation with reduced ordinary 
indices i', j' G {1, . . . , ni=i i^k n i} uniquely, that is it is the same as if we used appropriate 
versions of formulas (|T4|) . (|T5l) . 

Let us analyze, as an example, the case of a submatrix A' of matrix A of size N x N, 
N = nj=i n i- Let A' be built from rows of A indexed by % G A, with the order of rows 
preserved as is the habit when taking submatrices. Because for i,j indexing rows of A 
used to build A and for i',f corresponding to reduced vector indices for i and j there 
holds: 

i < j [i 1 ..i k - 1 ,a,i k+1 ..i r ] < [ji.-jfc-i, a,j k+1 ..j r ] <^=^> (19) 

[i 1 ..i k -i,i k+1 ..i r ] < [ji..j'fc-i, jk+i-jr] i' < f 

where G {1, . . . , XYi=i^ k n i) 

it is justified to say that the m-th row of A forming A' is the m'-th row of A'. Here m! is 
obtained from m in the composition of relations: m — > [mi..m k _i,a,rn k+ i..m r ] — > 
[mi..rn k -i 1 mk+i..m r ] — )• m' . Equivalently we could say that the ^..mn, a, mfc+i..m r ]- 
th row of matrix A is the [mi..rrik-i, rrik+i. .m r ]-th row of its submatrix A. 

Now we return to the main course of this lecture. From now on we silently assume that 
all Kronecker factors in (fT21 ) are of size at least 2. This we do for convenience of thinking, 
though we do not exclude the possibility that some or even all of the below arguments 
remain correct without this assumption. We will address this issue later after announcing 
the main results at the end of this and at the beginning of the next section. 

Using the introduced vector indices, we will split set M.u defined in (jTOl) into disjoint 
subsets {M v ) {hlMt .. tkp)t{vi>Va> .. >VpV 
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Definition 2.1 

,(Vl ,V2,..,Vp) 

(20) 

iki = hi = v x e {l..n kl } 
h P = jk p = v r e {i-n kp } 

where we assume that 1 < k\ < . . . < k p < r, 1 <p < r — 1, and 

(^)() )0 = 4^ = h^Ji, it^h, , (21) 

where () means the empty sequence of indices. 

The corresponding subspaces will be denoted according to the following 
Definition 2.2 

span ((Mu) Qi0 ) (22) 
span ((Ma) (fcli .. jfep)>1 ,.., Up) ) (23) 



[Mu, 



def 
'(),() = 

'(fci,..,fc p ),(i;i,..,ijp) 



Note that in both of the above definitions the notation refers to the assumed Kronecker 
product structure of U . Later we will meet U deprived of one or more of its Kronecker 
factors. Accordingly, the resulting U' will be assumed to be built from a fewer number of 
factors, and this will be reflected in the notation used for subsets of M.yi and subspaces 
of M.xji. 

Any vector in M.jj (see (Q) can be written as a linear combination of components each 
of which belongs to one of the above defined subspaces, hence Mu is the algebraic sum 
of all of them, 

M v = (24) 

(Mu) 0:{) + J2ke{l..r} J2ve{l..n k } ( M u)( k ),(v) + 
+ Z/fci<fc 2 e{l..r} £->(vi,v 2 )e{l..n kl }x{l..n k2 } ( M t / )(fci,fc 2 ),K,-«2) + ■■■ 
• • • + X^ 1 <..<fc r ._ 1 e{l.. r } /2(vi,..,v r -i)e{l..n kl }x..x{l..n kr _ 1 } (^)(fci,..,fc r -i),(j;i,..,-iv--i) 



Our aim is to show that ([24]) is in fact a direct sum, which will allow to estimate its 
dimension. In the meantime a little preparation is needed. 

Using notation similar to that of MATLAB, let, with y at the Z-th position, 

.,i r \,[n 1 ,..,n l _ 1 ,y,ni +1 ,..,n r - 1 ,n r 
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be a subrow of the [ix, ■■, i r ]-th row of matrix X of the size identical with the size of the 
Kronecker product (fl2l . The subrow is composed of entries of .., j r ]-th columns for 
which ji = y. The horizontal order of entries is preserved. 



Lemma 2.3 Let [ii,..,i r ], [ji,--,j r ] be vector indices of the Kronecker product KW\) cor- 
responding to ordinary indices i,j. 



a) If it ^ jk then any [bi, .., b r ]-th row of Ay , Sy satifies: 

E \ A u J) ] = 

— 1 L J [&i,..,6 r ],[:,..,:,c fc ,:,..,:] 



Cfe = l 

Eh 

c fc =l 



(ij) 
(7 



[6i,..,b r ],[:,..,:,c fe ,:,..,:] 







where Ck is at the k-th position. 
b) If i k = j k then 



E[ 

Cfe=l 



.4 



(ij) 



[fel,..,6 r ],[:,..,:,c fe ,:,..,:] 



[6i,..,6 r ],[:,..,:,Cfc,:,..,:] 



(i'j') 



c(i'j") 



[6i,..,6 fe _i,fe fe+ i,..,b r ],: 



[6i,..,6fc_i,6fe +1 , ..,6r],: 



(26) 
(27) 



(28) 
(29) 



for any [b±, .., b r ]-th row such that bk = ik = jk, Ck as above. On the right hand 
sides of (TJ2P and (TJP)) stand [£>i, .., bk-i, b k+ i, .., rows of matrices Ay,'^ ' , Sy,'^ ^ 



constructed from the Kronecker product KW\) deprived of its k-th factor: 

U' = [/« g> . . . ® C/^" 1 ) g) f/( fe+1 ) ® . . . ® C/W . (30) 

Ordinary indices %' , j' correspond to the accordingly reduced vector indices 
[h, --,ik-i, ik+u ••) v], [ji, ••, jfc-i, Jfc+i, • •,>]■ 



Proof 

a) is of course true for any row of , S^'^ indexed by a number different from i, j, 
which is in fact a zero row. As the i-th and j-th rows of Ay, Sy are of opposite signs, 
we show a) only for the i-th row. We do it by proving that every element of the left hand 
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side of ( 126|) and (|27|) is equal to 0. Indeed for My defined in ([6]) we have 



iik 



Cfe = l 



M, 



u 



[ll,..,l r J,[Cl,..,C r J 



(31) 



b'l >••>>], [ci,..,Cr] 

E?=i (1) k, ci • • • • • [^ (r) k J • (ic^U • • • • • P^U 

+ 0i , 

so its real and imaginary parts satisfy (1261) and (I27j) respectively. 

Now part b). If [fo l3 .., 6 r ] satisfying 6 fc = i fc = j k is different from [ii, .., i r ], .., j r ], then 
the ordinary index b' corresponding to the reduced vector index [bi, .., b k _i, bk+i, ••, b r ] is 
different from i', j'. In this situation both sides of (|28|) . ( 129]) are zero rows. 

Next we consider the i-th row, [bi, ..,b r ] = ..,i r ]. Again, we investigate the sums 
and (1291) entrywise. As above, we have: 

i>k 



Cfc=l 



(32) 



|*l,..,l r J,|Cl,..,Cr 

2 



+ ... + 



(nu ^ ^ (0 k C! ) • (nu ^^ w u,c 



hnr-> i fc-l> i fc+lr->MJ c lr->Cfc_i,C fe+1 ,..,C r ] 1 [ji ,. . ,j fc _ i J fc + 1 ,. . jV] , [ci ,. . ,c fc _ 1 ,c fc + 1 ,..,c r ] 

where the ordinary indices z', f, c' correspond to the vectors indices .., ifc_i, ik+i, ••, v], 
[ji, -Jk-iJk+i, -Jr], [ci,..,c fc _i,cj fe+ i,..,c r ] respectively. Now, by taking the real and 
imaginary part we get equalities (I28p and (T2"9"l) . Calculations for the j-th row are com- 
pletely analogous with the sign changed. 



The properties stated in the above lemma extend to linear combinations of matrices 
Lemma 2.4 Let B e (Mu)^^^. Then: 

a) For any k s E {hi, .., k p }, any Ck s E {l..nk a }, and any [b±, .., b r ]-th row in B: 



n k s 



[6i,..,6 r ],[:,..,:,Cfc.,:,..,:] 



\. u(t XJ ■ £ 



(33) 



d k .=l 



where c ka , d kg both at k s -th position designate subrows of row [B\< bi b i ,. 

In other words: the Ck a -th subrow is a multiple of the sum of all subrows designated 
by all the values of the k s -th subindex. 

b) For any k' £ {ki..k p } and any b-th row in B: 

n k> 

[B] b j. , c / , j = where c k i is at k'-th position. (34) 

c k' =1 



Proof 



The part a) of the lemma is of course true for all rows indexed by [b± , . . , b r ] such that 
b kl 7^ V\ or b k2 7^ t'2 or ... or bk v ^ v p , since these are zero rows coming from zero rows in 
4 J) , spanning (M^)^..^^..^. 

Let us thus consider the [bi, .., 6 r ]-th row of B with b kl = v\, b k2 = i>2, b kp = v p . Since 
B E ( M u) {hl ,..,k p ),(v u ..,v p ) for the [h, ..,b r }-th row we have, 



\B] 



[bi,-,b r 



(35) 



V (a) a&ti 



.4 



(id) 



s, 



[bi,..,b r ],- 
[6i,..,Mv 



+ 



and analogously for the Cfc s -th subrow designated by the value of the k s -th index, where 
we ommit the ranges for a^'^ and a^': 



l B }[b u ..,b r U-;-,-;C ks ,-;-,--] 



E....« 



(hi) 



A 



(iJ) 
U 



J [b 1 ,..,b r ],[:,..,:,c ks 

Hence, on the one hand: 



(i,j) 



(i,j) 

u 



\b 1 ,..,br\,[.,.. ! :,c ks , :,..,:] 



(36) 



,6 r ],[:,..,:,dfc ,:,..,:] 



e.... e;;:i , 



,4 



(7 



E....« 



(<j) 



[&l,..,iy]> [:,..,:, d ka ,:,..,:] 
[bi,-,b ks - 1> b hs+1 ,..,b r ] > : 



+ E....^ u) E£=i 



(id) 
u 



[bi,..,br],[.,..,:,dk B >-> 



- [6ii -,6^-1,6^+1, -,6f] 



where we used Lemma [2.31 b). which can be employed since the indices i,j in the sums 
satisfy i kl = j kl = vi, i kp = j kp = v p , and in particular i ks = j ks , as well as b ks = 
v s = ik s = 3k a - As previously, the ordinary indices i', j' correspond to vector indices 
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[i 1 ,..,ijfc,_ 1 ,i fc , + i,...,i r ], \ji,-,jk s -iJk s +i,-,jr) and designate nonzero rows in matrices 



Sjj,'^ constructed from the reduced Kronecker product 
Tj( k ° +l) ® .. ® U {r) . 

On the other hand the expression ( 156)) translates into 



[/■(*.-!) 



71 C/' 



> c fc s 



(i'j') 



Djj, 



[b 1 ,..,b ks _ 1 ,b ks+1 ,..,b r ],: 
[6i,..,6fc s _i,6j, s+ i,..,6 r ], 



+ 



with [/', i', j' described above. Combining ( 157)) and (J5HJ) we obtain a). 

To prove the part b) we use again (I35p and rewrite it for the subrows corresponding to 

k' t {h..k p y. 



b r l, [:,..,:, c,./,:,..,:l 



.4 



(7 



J [&i,..M 



5. 



J [6i,..,fc r ],[ : 



(39) 



Here i, j are such that iy ^ jy- We may thus use Lemma [2.31 a) to find that summation 
of subrows (1591) produces a zero vector: 



Yl [■ B ][6i,..,6 r ],[:,..,:,« V ,:,..,:: 



c,./=l 



E....« ( ^ES'=i 



.4 



[6l,..,6 r ],[:,..v,c fc ,,:,..,:] 

V ■ + V trW) -0 = 



5, 



[6l,..,6 r ],[:,..,:,c fc ,,:,..,:] 



Related is the following result which we will also need to prove the theorem on a direct 
sum that comes next. 

Lemma 2.5 Let [Mu] be the space spanned by matrices \ S^' such that = 

V / k,v 

jk = v, where v E {l-.n^}. There exists an isomorphism mapping [Mu) onto Mjy/, 

V / k,v 

where U' = <g> . . . <g> ® <g> . . . <g> t/ (r) . It can be chosen in such a way that 



any subspace (Mp)j tl ^ v ^ of yMu J wi/i k q = k and v q = v for some q £ {l..p}, 



zs mapped onto the subspace (Mu')u^ ..^A, (vi, ..,v q -i,v q+ i, ..v p ) °f M w , where k[, k' v _ x 
are positions ofk\, fcp m i/ie sequence (1,2, ...,k — 1, fc + 1, r). fiVoie 

t/iat fcg indicate positions of factors U^ kl > , U^ kq - 1 \ U ( - kq+1 \ in the product 



U' .) In particular the chosen isomorphism maps (Mj/)/^ ^ 
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onto (M[//' 



0.0 ' 



Proof 

We will show that matrices sjj'^ spanning (^Muj remain in a simple relation 

with matrices Ay/* \ S^,'-* \ Here i', j' correspond to vector indices .., ik-i, ifc+i, ••, i r ]> 
b'i) ••) iit-ijife+i) ••) jr] °f the reduced Kronecker product [/' (I30I) . i.e., are given by the 
appropriate versions of (fT4l) and (TT5|) . 

We will be more precise about the above mentioned relation later, but at first let us state 
it without too many details. Any spanning [Mu] contains, as its only nonzero 

and disjoint submatrices, copies of the corresponding A^,'-' ^ , where each copy is 
multiplied by one of rik multipliers. Column positions of the copies do not depend on 
The multipliers, of which at least one is nonzero, depend only on the column positions 



of the copies they act on. The same applies to and S^,' J 1 , with the pattern of 

submatrices and values of the multipliers identical with that for A^} and A^ \ 

For a more detailed explanation let us consider the z-th row of A$ j \ The 

[di, .., dk-i, dk+i, .., rf r ]-th element of the dk-th subrow of it, dk G {l..rik}, reads: 



which is equal to 



.4 



i,[di,..,d k ,..,d r ] 



Re 



V 



i,[di,..,d k ,..,d r ] 



Re ([^]i 1 [di,..,d fc ,..,d r ] ' Wlj^du..,^,..,^]} 



(41) 



Re ([UiX 



di 



31, di 



.... [*/<*>] 



ik.,dk 



3k, dk 



.... [U^] 



where lU^] . , 
a result we get 



- 3k,dk 



can also be written as [f/w] v d 



v ,dk 



jr jC^ 

[[/(*)] 



, (42) 



v,d k 



. As 



■u.dfc 



is A 



2 / 

]v,dj. ' [J^ \h,--,ik-i,ih+l,--,irUdi,--,dk-i,dk+i,--,dr] ' [j u ..,j k _ 1 ,j k+1 ,..,j r ],[d 1 ,..,d k _ 1 ,d k+l ,..,d r 



Re 



i',[di,..,d fe _ 



l:"*:+l!--! a r) 



u,d fe 



A (i ' J,) 



i',[d 



l.-«-l«+lr 



The equality of the far left hand side of ( l4Tj) and the far right hand side of ( )43|) is an entry 
by entry equality between two rows. Now, if we change the sign of Re in the expressions 
in (JHJ, (142|) and (143|) we will obtain a relation between the dk-th subrow of the j-th row 
in A^ ,j) and the j'-th row of Ay,^ \ If we used Im instead of Re, we would get a relation 
between the dk-th subrow of the z-th or j-th row of and, correspondingly, the z'-th 
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or j'-th row of S\)^ \ This is summarized below: 



A (id) 






= A*fe • 


71 t/ 


AiJ) 

A u 


j,b, 


■vAr.-v 


= Aifc ■ 




q(id) 






= Hk ■ 

] 


'q{i 


q(id) 




..,:,dfc,:,..,: 


= /^fc • 


'q(i' 



(44) 



where fi k 



[U^] vd ■ One of the n k numbers \x k is nonzero since is unitary. 

to denote a submatrix of X built of those x-th rows 



Now let us use [^] [ . i .. i . i6fc) . j .. ) . ]>[ . j .. ) . ;Cfer 
and y-th columns of X for which x k = b k and y k = c k . As a consequence of relations 

(144j) we see that , spanning (Mt/j are all zero matrices except for their 

submatrices: 

,:,v,:,..,:],[:,..,:,c k ,:,..,:] 



A 



(id) 
U 



(id) 

u 



where /ik 



c k e {1, ..,n k } , 
c k E {1, ..,n fc } , 



(45) 
(46) 



(where w is at k-th position), with one of \i k being nonzero. Analogous equalities we 
could write for linear combinations on both sides. ( (I45j) and (l4*6l) can be checked entry- 
wise, in which one uses that the [bi, .., bj~-i, bk+i, ■■, b r ] , [c\, .., Ck-i, Ck+i, ••, c r ]-th element 

:,6fe,:, ..,:], [:,..,:, c fe ,:,..,:] * s ec L ua l to [-^][6i,..,6 r ],[ci,..,c r ] ' Here of course we identify, for ex 
ample, [bi, .., bk-i, bk+i, ■■, b r ] with ordinary index V . See also the comment that follows 
the introduction of vector indices.) 

Since all \i k add up to 1 , i.e. the norm of the v-th. row of U^ k \ we can choose an 
isomorphism \l/ between ( M v j and M.u> to be defined as (again assuming that v 

V / k,v 

stands at fc-th position): 
so that 



k,v 



E 

c fc =l 



Ai • 



(47) 



VIA 



VIS; 



4 li 



(*',/) 
17' 

(i'J') 



(48) 



Consider now (Mu)^ kl kp ) ( V1 v p ) c (-^7) with k q = k and v q = v for some g G 



{l,..,p}. Matrices S$ J) spanning this space are indexed by z,j for which i kl 



k.v 
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j kl = v u i kq = j kq = v, i kp = j kp = v p , and finally i, ^ j, for / ^ {k 1} .., k p }. Note 
that according to (148]) the isomorphism \1/ maps these matrices into all matrices spanning 
(Mr//)/,, \ / w where .... fc' i are defined in the lemma. This is 

because if we take any pair z, j designating , 5"^ spanning (Mrj)( fel kp ^ ^ „ ), in 
the reduced vector index .., ifc_i, ifc+i, ••, i r ] a t positions fc^, k' v _ x we have subindices 
ifci, ifc,+i, ••• ifcp that are equal in pairs with subindices j fel , j kq _^ j kq+1 , 

... j kp at the same positions in the vector index [ji, .., j k -i, j k +i, .., j r }- Their values are 
V\, v q _i, v q+ i, v p , respectively. The only other pair of equal subindices, namely 
ik q — jk q — v, is absent from the reduced vector indices corresponding to At other 
positions we get any possible pair of different values as we take any pair i,j allowed here, 
that is indexing A$ 3 \ S^ 3) spanning (Mu) iki ^ kpUvi ^ Vp) . ■ 

Our most important result in this section is the following 

Theorem 2.6 The algebraic sum of subspaces in (Efy is a direct sum. 

Proof 

We will prove the following statement equivalent to the above theorem. 
Let X {)>() G (Mu) (U) , X (fel ,.., M 

,(Vl,..,Vp) ^ (Mrj) 

(fci,..,fcp),(«i,..,i;p) ^ or a ^ possible 
choices of (fci, .., k p ) and (v±, .., v p ) be matrices whose sum is equal to a zero 
matrix. Then all these matrices (we will call them X-matrices) are zero 
matrices. 

The proof will be by induction with respect to r, the number of factors in Kronecker 
product (PD - 

If we have only one factor, U = t/ (1) , then (Mrj) () () = Mrj is a direct sum of one 
component, namely (M[/)q q. 

Assume now that the theorem is true for any (r — 1) factor Kronecker product (|T2l) . In 
particular that what we stated above as eguivalent to our theorem is true for such a 
product. 

Let our X-matrices add up to a zero matrix and let one of the summands, Xt kly _ jk ) 1 ( vi ,..,v ) 
be nonzero. Let -Xa liiii y i ( Wri , ) be simply denoted by X. Thus X is the sum of the 
remaining X-matrices. We also need to introduce k[ < k' 2 < ... < k' r _ p such that 

{k l ,..,kp}U{k' v ..X- P } = r}. 

The remaining X-matrices, forming X in a sum, can be split into two groups: 

• Into one group we put these belonging to subspaces (Mrj)^ l ^ Wq) such that 
k[ {li, .., l q }. ( An X-matrix belonging to (M[/)q a is in this group, if {ki, .., k p } 
is nonempty. The other situation we consider at the end of the proof.) Thanks to 
that choice we can apply Lemma 12.41 b) to all of them, as well as to their sum, 
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which will be denoted by Y: 



Em 



for any 6-th row 



(49) 



Into the second group we put all the X-matrices which belong to (My)^ y / 1 ^ 
such that k[ G {l±, .., l q }, that is l s = k[ for some s G {1, .., q}, and then w s — v G 
{1, .., ny }. The group can further be split into subgroups corresponding to different 
values of v. When we add all the matrices in the whole group to form their sum 
Z, the 6-th rows of Z, where the corresponding [b\, ..,b r ] satisfy = v, are sums 
of 6-th rows of matrices belonging only to the subgroup associated with value v. 
Matrices from other subgroups have at these positions zero rows. 

We can apply Lemma [2.41 a) to each matrix in a subgroup, as well as to their sum 



[Zv] 



(50) 



Since, as we have just said, matrices forming Z v in a sum are responsible for for- 
mation of only those 6-th rows for which by = v , we can write for the sum Z of all 
Z v : 



[Z] 



Mr 



K'. fih 



\z\ 



(51) 



At this moment recall that, using our symbols, X = Y + Z . Also note that if one had to, 
one would classify X as belonging to the first group since X G (Mu)( kl fe ) ( Vl Vp ^ and 
k[ {ki, .., k p }. So X satisfies (l4"9~j) with Y replaced by X and analogous property holds 
for Z through the linear combination Z = X — Y . The last statement implies that the 
sums on the left hand side of fl5Tj) are zero rows, which in effect gives Z = 0. Further, 
let us return to component Z v of Z. As Z is a zero matrix and Z v have nonzero rows in 
disjoint regions, all Z v must be zero matrices. 

Recall that Z v are formed by X-matrices belonging to all possible (Mjj)^ z ) ( Wl w ^ 
being subspaces of the larger space ( Mu ] introduced in Lemma 12.51 and spanned 

V / k[,v 

by A$ j \ S$ j) with 

iu = 3k' = v - Lemma [23] tell us that this larger space is iso- 
morphic to M w where U' = C/W <g> U^'^ <g> <g> C/W. We denote this iso- 
morphism by that is ^ ^My^ ^ J = My/. In accordance with Lemma 12.51 we 

choose isomorphism \1/ such that it maps any space (M^)^ ^ ^ ^ under consider- 
ation onto space (Mm)/,, „ \ , x where, as above, L = k\, w s = v, s G 
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{l,..,q}, and where l[, V x are positions of l\ Z s _i, l s+ i, l q in the sequence 
(1,2, ...,k[ -l,k[ + l,...,r). In particular, (Mu) {ls) ^ Ws) = (Mu)r k ,\^ v) is mapped onto 

(M^)(), - 

The images under \l/ of X-matrices forming Z v are contained in the above mentioned 
image spaces (Mjy/)^, v \ ^ ^ itu +i which in turn are contained in M^/. The 

images add up to a zero matrix since the very X-matrices forming Z v add up to a zero 
matrix. 

On the other hand, and this is our induction assumption, M^/ is a direct sum of 
(My/)^, ; , \ w lW+l Wq y for we have all possible parameterizing pairs of sequences 

coming from all possible pairs (li,..,l q ), (wi,..,w q ) such that l s = k[, w s = v for some 
s G {!,..,<?}. 



( According to Def. I2.1[ sequences (li,..,l q ) as well as are ordered. Let se- 

quences (liT-,l q -i) and (wx, .., w s+ i, .., w q ) be given, with entries in appropriate 
ranges. Here s is defined to be such that l[ < . . . < < k[. Then for ..,l q ) = 

l' s _ lr k[, l' s + 1, . . . , 1'^ + 1) and [w x , . . . , w s _ u v, w s+1 , .., w q ) space {M v ) {lu „ :lqUwi> „ >Wq) 
is mapped onto (Mn-/)/,, „ \ , v ) 

Thus all the images of X-matrices forming Z v are zero matrices and the same must hold 
for the very X matrices. As v is taken to be an arbitrary element from the set {1, .., ny }, 
all matrices forming Z, that is belonging to the second group, are zero matrices. So X 
is formed only by matrices from the first group. 

Analogous reasoning can be repeated with k' 2 , k' r _ p instead of k[. As a result we 
obtain that X is formed only by X-matrices belonging to (Mj/)^ j j , Wi Wq) such that 
k[, k' r _ p do not belong to {h,-.,l q } or, in other words, by X-matrices belonging to 
( M u)( ko ^.., kgt )A W i,,w t )i where { k 9i>~i k 9t} C {fci, .., 



If one of these , call it X', is nonzero and belongs to (M[/)/ fe fc w Wi with t < p, that 

is .., k gt } 7^ {/ci, .., /c p }, it can be expressed with the use of the remaining ones and X. 
We can repeat the whole reasoning with X' playing the role of X and k[ being replaced 
by some k! G {1, .., r} \ {k gi , .., k gt } and simultaneously k! G {k±, .., k p }. We then classify 
summand X as belonging to the second group of matrices, since X G (Mu)( kl fc ) ( V1 „ ) 
and fc' G {/ci, .., fc p }. We show as above that this second group of matrices consists of 
zero matrices which contradicts our assumption that 1^0. 

Note that at this stage (Mu)q q could appear as (Mu)^ fcg ^ ^ wi m y The short version 

of our reasoning for X (the role of X' at this stage) belonging to (Mu)q q is provided at 
the end of the proof. Also, precise repetition of our reasoning is always possible because 
no components in matrix sums are missing - even if we write that a matrix is formed only 
by some other matrices this means that the remaining ones are there as zero matrices. 

Having X G (M[/)/ fci fc > , ^ on the left, on the right we can only have nonzero 
components belonging to (Mu)^^^^ with (w x ,..,w p ) ^ (v^.^Vp). Obviously 
(wi, ..,w p ) are also distinct, since any component (Mu)( kl k ) ( m w ) in what we aim to 
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prove to be a direct sum, occurs only once. Collecting both sides together we get that the 
sum of X-matrices (—X is one of them) belonging to (Mu)( kl k > ( Wl w y for sequence 
(ki, .., k p ) fixed by the choice of X and for all possible sequences of values (wi, .., w p ), is 
a zero matrix. Since regions where these last candidates for nonzero X-matrices have 
nonzero rows do not overlap, also they must be zero matrices. This is exactly what we 
wanted to prove, hence indeed ( 1241) is indeed the direct sum. 

Now, as promised twice in this proof, we consider the case of X £ (M[/)q q, that is the 
case of {ki, ..k p } being empty. Splitting the remaining summands into two groups, as 
previously, for any k' £ {k[,..,k' r } = {1, ..,r}, we conclude that any X-matrix belonging 
to (My)^ i)( Wl w ) such that k' £ {li,..,l q } is a zero matrix. Taking all k! under 
consideration we find that all X-matrices forming X must be zero matrices. So X must 
be a zero matrix too. 



Let us rewrite (I24p in a direct sum fashion: 
M v = (52) 

(McOo.0 © © fc6{1 ..r}©«6{l..» fc }(^) (fc) , (w) © 

© ©*i<fc 3 e{i..r}©(t» 1 ,« 3 )e{i..t» fcl }x{i..»» a } ( M u)(ki,k 2 ),(vi,v2) ® ••• 

••• © ©fci<..<fc r -ie{l..r} ®(vi,..,v r -i)e{l..n kl }x..x{l..n kr _ 1 } ( M ^)(fci,..,fc r _i), (vi,.. ,v r -i) ' 

Calculation of the defect of U requires calculation of the dimension of the Mjy, which 
luckily can be decomposed into the dimensions of the direct sum components. We will 
make this decomposition more informative in a moment. 

3 Estimating the defect 

Calculation of all the dimensions of the direct sum components in ( 152]) can be reduced 
to calculation of dimensions of subspaces like (Mu)q q considered for subproducts of the 
Kronecker product U in ( IT2l . 

Lemma 3.1 There holds: 

dim ((MtOcfa,..,^,..,^) = dim ((M^p) () , () ) (53) 

where 

IP = (g) U {k) (54) 

fe6{l..r}\{fci,..,fep} 

zs the original Kronecker product (TT^j deprived of its k\-th, k p -th factors. 
Proof 

In the proof we use Lemma [2.51 repeatedly, assigning appropriate values to U, k and v. 
We have the following relations: 
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Let U = U, the initial Kronecker product f[T2"j) . and let k = k p , v = v p . Then 
by Lemma[221(M [ /) (fei) )fcp) (ui ) is isomorphic to (MaO^,..,^),^,..,^), where 
£7' = 0; e {i r }\{fc \ U">. In this case k' p _ 1 mentioned in Lemma 12751 have 

values k[ = k±, k' 1 = fc p _i, as they are the positions of ki, k p _\ in the 
sequence (1, 2, k p — 1, k p + 1, r). 

Let U = U' , k = kp-i and u = Then (Mj//)^ fc ^ v _ > is isomorphic to 
(Mt/")^,..,^),^,..,^), where £7" = (8>^ {1 .. r}X{fep _ 1)fep} V«. 



• Let £7 = f^ 1 = ® le{1 „ r} \ {k2 „ kp} U®, and let fc = h and v = v v Then (M^-O^),^) 
is isomorphic to (Myj^Q, where U ? = (8)ie{i..r}\{fci..fcj,} 

Thus both subspaces mentioned in the lemma are isomorphic through the composition 
of the above listed isomorphisms, so they are of equal dimension. ■ 

Let us simplify our notation by introducing, for U being the product (|12|) : 

= f dim ((M^) (),()) hi U> = (g) 17« (55) 

fce{l..r}\{fci..fcp} 

4 d ^ f dim ((M [/ ) ())() ) (56) 

Using the fact that Mjj is the direct sum (I52"j) and that the dimension of any component 
space (Mu)r kl k \ r vi v \ can be calculated by Lemma [3TT1 to be d^ ki kp y, the dimension 
of M[/ can be expressed as in the Theorem 13.21 below. 

At this point note also, that for a (r — l)-element set {ki, k r _{\ = {l..r} \ {s} there 

holds <i| / fci kr ^ = dim ^(M^c^q qJ = dim (Mj/w). We do not assume any Kronecker 

product structure for CA% that is vector indices are identical with ordinary indices here. 
Using such vector indices we cannot fix their only component, (-Mu( s ))(i) ( v ) would be 
empty. 



Theorem 3.2 Let U be a Kronecker product of unitary matrices (TB) . Then the dimen- 
sion of the space fAjj constructed for U according to (TJ|) reads: 

dim (Mu) = (57) 

d{} + J2ki€{l-r} n fei ' ^{fei} + ^2k 1 <k 2 e{l..r} n ki n k 2 ' dfkxM} + 
+ Sfc 1 <...</c r -iG{l..r} n k 1 ■ ■ ■ ^fc r _i ' d{ klt kr l y 

Accordingly (see Definition ll.ip . the defect of U equals to (77 — l) 2 — dim (Mf/). 
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Maybe it is a good moment to wonder about applicability of our formulas when there 
are lxl Kronecker factors in product (fl2|) forming U. This has been promised in the 
paragraph preceding Definition 12.11 

First, recalling this definition, consider subset (M-u)^ k 2 k ) Oi v 2 v p ) °f se ^ -M-u- If k G 
{1, ..,r}\{k\, ..,k p } and n k = 1, this subset is empty The reason is that you cannot find 
any Ay , belonging to it, that is indexed by i,j satisfying i k ^ jk, because i k , jk G 
{1, ..,n k } = {1} . In this case let the corresponding subspace (My) ^ 2 k ^ v ^ v ^ of 
space Mf/ be defined as the zero space {0}, which is a space of dimension 0. In particular, 
when there are lxl factors in product ffT2|) . (M.u)q a is empty and (Mu)q a we define 
as the zero space. Consequently, d^ kl kp y has value 0, if for some k G {1, .., r} \ {fci, .., k p } 
there holds n k — 1. In other words U', used in the definition ( )55|) of d^ ki kp y, must not 
contain any lxl Kronecker factors, if this d^ ki k y is to be considered as a potentially 
nonzero number. (It still can be zero if for example U' is a permutation matrix, in which 
case all A^'r, S^, are zero matrices, so (Mu')q q is the zero space.) 

Now, let us remove in sum (I57j) all components d^ ki fc y that must be zero due to the 
presence of 1 x 1 factors. What is left are only those d^ ki k y, for which there are no l's left 
in sequence (ni, ..,n r ) after removing the fci-th, k p -th entries. The sum ( 15T|) truncated 
in this way looks as if it were written for U deprived of its 1 x 1 Kronecker factors, in 
which case we accept the validity of formula f )57|) . This new U differs from the original 
one by a unimodular factor, so the space and the defect d(U) remain the same, as 
one easily finds noting that first of all A^' 3 , S^^ do not change with multiplication of 
U by a unimodular number. The answer to our question is thus: yes, we can use (|57|) in 
presence of 1 x 1 factors, provided that appropriate d^ kl kp y are assigned O's. 

This has some practical importance. For example, if we wanted to calculate defects for a 
large number of Kronecker products of at most r factors, we could store the sequences of 
sizes as rows of a matrix with r columns, filling rows with a certain number of l's where 
the number of factors would be less than r. These rows would then be the input data for 
a procedure calculating the defect of an r factor Kronecker product of unitaries. 

We can go even furher. Let us also, in this context, return to the direct sum formula 
(1521) . We have adopted above a convention, according to which (Mi/)r kl k2 k \ ( V1 V2 v ) 
is the zero space if there are l's in sequence (rii, .., n r ) after removing from it the ki- 
th, kp-th entries. Let us throw out these zero spaces from fl52|) . We are left with 
(fci,fc 2 ,..,fcp), (vi,v 2r -,v p ) wriere an k associated with n k — 1 sit in {ki, .., k p }. But every 
such (Mr,)^^^.^) can be replaced by an equal space ^ud^^^^^ > 
where: 

• U' is U deprived of its 1 x 1 Kronecker factors, with the order of the remaining 
factors preserved, 

• (k[, k' 2 , .., k' pl ) contains the positions, of subsequent entries of the subsequence ob- 
tained from (ki, .., k p ) by throwing out all fc's such that n k = 1, in the subsequence 
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obtained from (1, 2, .., r) in the same way (This is cumbersome, being a consequence 
of our notation, and we met such formulation in Lemma [2.51 ). 

• (v[,v' 2 , ..,v pf ) is the result of throwing out the entries (of the only possible value 1) 
out of (fi,f2, ■■jVp) corresponding to entries k of (k±, .., k p ) for which nk = 1. 

(Both spaces are equal because they are spanned by the same matrices. In general, 
= A^; 3 \ where i',f correspond to the reduced vector indices obtained from the 
vector indices corresponding to i,j, respectively, by removing those fc'th entries of these 
vector indices (entries necessarily equal to 1), for which k satisfies rik = 1.) 

After this truncation and replacement ( 1521) looks as if it were written for U', in which 
case we accept the validity of the direct sum formula because U' does not contain lxl 
Kronecker factors. Since Mu = My, the only difference between U and U' being a 
unimodular factor, the starting point direct sum with all those zero spaces appears to 
lead to the correct direct sum we would write if we first removed the lxl Kronecker 
factors. So, we can go either way, removing or leaving these factors. 

It is interesting to give the lower bound for the defect of U being a Kronecker product, 
because, as it was pointed in the Introduction, for most unitaries the defect is zero. The 
lower bound we give here is associated with upper bounds on the dimensions d^ ki k ^ of 
(M{//)q q for the Kronecker subproducts U' as defined in (I55I) . all of which are in fact 
bounds on d^ with U replaced by an appropriate U' . All we need to know about these 
bounds is contained in the next lemma. Note that although we assume that the sizes of 
Kronecker factors are greater than 1, any factor of size 1 - a unimodular number - can 
be absorbed into one of the factors of size greater than 1. 

Lemma 3.3 Let U of size N x N be the Kronecker product [W(l with r > 1 and factors 
U( r ' of size ni x n\, n r x n r , where all nk > 1. 

Then the dimension of the space (Mu)q q constructed for U (Definition ^. tfy is bounded 
in one of the following ways: 

• If nk > 2 for any k G {l..r} ; then 

4 < (iV-l)(m-l)-...-(n r -l) . (58) 

• If = . . . = nk p = 2 for some distinct values k s £ {l..r} ; and for k (jt {k\..k p } 
there holds > 2, then 

d{ } < (N - 2 p ~ 1 )(n 1 - 1) • . . . ■ (n r - 1) . (59) 

Proof 

From Lemma 12.41 b), if B £ (M{/)q q the subrows of any of its rows satisfy, for any 
k £ {l--^} (where is at fc-th position): 

E ^M:,..,^,,,] = > ( 6 °) 

Cfc=l 
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as in this case the role of the set {k±, k p } from Lemma [2.41 is played by an empty set. 
We will show that the b-th row of B can be parametrized by at most {n\ — 1) • . . . • (n r — 1) 
independent parameters. For example, let [B] b , Ci i be known for c\ G {1, .., {ri\ — 1)}, 
c r G {1, .., (n r — 1)}. (J60l) determines the remaining entries in the following order: 

# [ B k[n u c 2 ,..,c r ] = ~ Ec^T 1 ' [ B k[c u c 2 ,..,c r ] t0 make a11 l B }b,[c u ..,c r ] kllOWn fOT C l ^ 

{1, c 2 G {1, .., (n 2 - 1)}, c r G {1, .., (n r - 1)}. 

# l B }b,[ Cl ,n 2 ,c 3 ,..,c r ] = -T}c 2 =l l) [ B }b,[ Cl ,c 2 ,..,c r ] t0 make a11 l B }b,[ Cl ,..,c r ] kllOWn fOT C l £ 

{1, c 2 G {1, ..,n 2 } c 3 G {1, .., (n 3 - 1)}, c r G {1, .., (n r - 1)}. 



• \B] h < „ „ „ „ i = — y^ nr i f-Bl, r „ „ „ , to have all the entries of \B], 

L J0,[ci,C2,C3,..,C,--l,nrJ Z— /Cr = l L J 6, [ci ,c 2 , . . ,c, — \,c r \ L J o,: 

known. 



Next we show that the 6-th row of B filled using the above algorithm, starting from those 
initially known (m — 1) ■ . . . ■ (n r — 1) entries, satisfies ( 16T)|) for any fc G {l.-r}. That is, that 
the following equality holds for any k G {l-.r} and c s G {l..n s } for any s G {l-.r} \ {A;}: 

\- B ]b,[c 1 ,..,c k _ 1 ,n k ,c k+l ,..,c r ] ~ ^ y [^]b,[ci,..,c fc _i,d fc ,c fe+ i,..,c r ] ' (^^) 

4=1 

Let in the above expression = n^, Cfc p = n fcp and c s ^ n s for s G {l-.r} \ .., k p } 
with ki < k 2 < ■ . . < k p . By this we mean that = Ck = c^ q = n kq for some q G {l.p}- 
(Do not mistake ki, fc p for those from the second item of the lemma.) Then the left 
hand side of (|61~]) . where some of the entries are obtained using the algorithm from the 
dotted list above, reads: 



^pK p -l) r Rl _ /_ v^Kp-i- 1 ) r R l 

^%=1 L jD Jfe,[ci,..A p ,..,Cr] ^<4 p = l ^ Z^d fep _ 1= l ^J{ ) ,[ Cl ,..A > _ 1 ,..A p ,..,Cr 



2^d kp =i y 2^d fepl =i ^•••^ 2^ 2 =i ^ 2^d kl =i l n lb,[c 1 ,..,dk 1 ,-,dk P 

i^d kq =i y 2^d kp =i y- y z^d kq+1 =i y 2^d kq _ l= i y- y 2-^=1 l-oj^ci,..,^,..,^ 

= - E2j=i [ 5 ] 6) [, 



I c\, ..,n kl , ..,n kq _ 1 ,..,d kq ,..,n kq+1 ,..,n kp ,..,c r ] l 



Z-A 9 =l I J6,[ci,..,c fel ,..,c fcg _ 1 ,..,d fe9) ..,c fcg+1 ,..,c fcj ,,..,c r ] 

En fc -1 rpl 
d fc =l L- D J6,[ci,..,d fe ,..,c r ] ' 

and equals thus to the right hand side of (I5TT) 

In this way each row of B G (M[/)q q could be potentially parametrized by no more 
than (m — 1) ■ . . . ■ (n r — 1) parameters - if there were no other restrictions caused by the 



21 



structure of - but we need to parameterize only N — l rows in general. This is caused 
by the property that all column sums in B are zeros, just as is the case with the spanning 
matrices A^\ Hence, in general, we have no more than (TV — l)-(n\ — 1)-. . .-{n r — 1) 

free parameters to determine B. In other words: all B £ (M{/)q q belong to a single space 
of dimension (N — l)(ni — 1) • . . . • (n r — 1). 

The case when some n k are equal to 2 is a separate one, but what really counts is when 
there are more than one n k equal to 2. 

In this case let again B £ (Mu)q q and let = 2, n kp = 2. The [&i, ..,6 r ]-th row 

of B is a linear combination of [&i, ..,b r ]-th rows of 4 <J) such that, amone; other 

conditions, i kl = 2/j kl , i kp = 2/j kp , where of course i kq , j kq £ {1,2}. 

Consider the group of rows in B indexed by such [b\ , . . , b r ] that for some chosen values 
si, s p £ {1, 2} the vector index satisfies 

(b kl = s 1 , . . . , b kp = s p ) or (b kl = 2/s 1 , ... , b kp = 2/s p ) . (63) 

Observe now that the rows of B from this group are linear combinations of respective 
nonzero rows of only those A^\ (from the collection of those spanning (M[/)q q) 

which are indexed by i,j satisfying: 

[i*i,»,ifcp] , [jki,-,jk P ] e [2/si, ..,2/sJ } . (64) 

Since subsets of (z, j)'s (of the set of all (i,j)'s indexing , spanning (Mu)q a) 
fulfilling condition ( |64|) for fixed si,...,s p are disjoint, also disjoint are the subsets of 
(<Mu)qq of those corresponding A^\ Sjj which are nonzero matrices. Within each 
such subset, associated with (|64l) . the nonzero rows of its members hit the index area 
defined in (163]) . Therefore the matrices A^ 3 , S 1 ^'^ belonging to it not only have the 
property that their rows add up to a zero row, but also that all the rows of such a matrix 
indexed by [bx, .., b r ] satisfying (153"]) add up to a zero row. 

Thus, through a linear combination, in the considered group of rows in B one of the rows 
depends on others. We have 2 p /2 = 2 V ~ X possible choices of [si, .., s p ], hence 2 V ~ X groups 
(of the specified type) of rows in B with 2 P ~ 1 dependent rows altogether. Using the 
previous method of parameterizing rows, we can parameterize no more than (TV — 2 P ~ 1 ) 
rows, introducing no more than (rii — 1) • . . . • (n r — 1) parameters into each row. In 
other words: B being contained in the described above (N — 2 p ~ 1 )(ni — 1) • . . . ■ (n r — 1) 
dimensional space is a necessary condition for B being a member of (Mt/) A q- 

Our final remark will concern the case when we deal with a one factor Kronecker product, 
r = 1, U = [A 1 ), N = ri\. The bound, either formula (]58j) or formula 05 9p . takes then 
the form (N — l) 2 . This agrees with these facts: 

• dim((M[/)Q q) = dim(M[/), no Kronecker product structure assumed for U . 

• Mj/ is contained in the space tangent to the set of all N x N doubly stochastic 
matrices, that is the space of all real N x N matrices in which entries in any row 
or column add up to 1. This space is of dimension (N — l) 2 . 
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In what follows Theorem 13.21 there are arguments behind defining (Me/-) a a as the zero 
space if there are lxl Kronecker factors in (1T2|) not absorbed into larger factors, and 
that dVy has to be assigned then. In this case we can also use the bound (1551) or (159]) 
which produces the correct dimension 0. 

Having the above result of Lemma 13.31 let us replace the values dV ki fc i. in Theorem 13.21 
by their upper bounds being the appropriate right hand sides of ( 155]) and (159]) . We will 
get an upper bound on dim(Mj/), equivalently a lower bound on the defect d(U) = 
(iV-l) 2 -dim(M c/ ). 

What is more, this can be done successfully also in the presence of unabsorbed lxl 
Kronecker factors. Since the bounds (l55j) . (1591) yield for d^ ki k ^ corresponding to zero 
spaces in the direct sum ( 152]) . whose presence is associated with the lxl factors, bounding 
the dimensions of these spaces does not affect the total upper bound on dim (Me/) at all, 
that is we obtain the value we would get if the lxl factors were first absorbed into larger 
factors. This can be better understood if one carefully follows the reasoning following 
Theorem 13. 2 1 The practical consequence is this: if one had some other bounds on d^ ki k 
yielding zeros for spaces that are necessarily zero spaces because of 1 x 1 factors, one 
could leave the lxl factors unabsorbed. 

To better expose properties of our bound, it is convenient to introduce another notion, 
the generalized defect D(c7) of a unitary U, defined as 

D([/) = d(U) + (2N-l) = iV 2 - dim (Me/) = dim ((Mf/) 1 ) (65) 

When we recall expressions in (j5J) and the preceding formulas it will be clear that D(t/) 
is the dimension of the space ([I]), that is the space: 

Vi/ = {R : \RoU = EU for some antihermitian E} (66) 

We use in our new definition the word generalized because the definition of d(U) is 
suited for the class of unitary matrices which have the dimension of manifold {D r ■ U ■ 
D c : D r ,D c unitary diagonal} equal to 2N — 1. Because of some applications of the 
defect mentioned in the Introduction it is more convenient subtract the dimension of this 
manifold from D(c7) to define the defect associated with other type of U. 
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The corresponding lower bound on T)(U) will be, thanks to Theorem 13.21 and Lemma I3TB1 
and definitions (|55|) . ( 156|) . : 



(ni • . . . • n r ) 2 - 

n „ t - 2«^ i -»-m . n (n,-i) 

^£{l..r} / ze{i..r} 

e in fel j n - 2^ 1 ->\W)- 1 ] • n (»«-!) 

fcie{i..r} \ ye{i..r}\{fei} / ie{i..r}\{fci} 

e n n, - 2^>\^.*>>m ] . n ( ni _i)] + ... + 

fci<fc 2 e{l..r} \ \/e{l..r}\{fci,fc 2 } / le{l-r}\{ki,k 2 } 

E K 1 -^- 1 ( II n, - 2^i..m{^-x})-i ) . "Q (n| _ i) ] ] ,(67) 
fci<...<fc r _iG{i..r} y ye{i..r}\{fei..fe r _i} J ie{i..r}\{fci..ftr-i} 

where for a nonempty subset ^4 = ...,/ s } of {1, ...,r} expression jj 2 ^4 denotes 

• the number of 2's in sequence (n^, ...,n/ a ), if there are any, 

• 1 otherwise. 

Although the above expression looks very complicated, below we show that it is, almost, 
a mere product of trivial functions of the sizes rii, n r of factors of the considered 
Kronecker product. 

Because in Lemma 13.31 it was assumed that all Kronecker factors are of size greater 
than 1, formally (16T|) is a valid lower bound only in this situation, therefore we confine 
ourselves in the formulation of the theorem below. But, in the paragraph preceding the 
one containing the definition (|65l) of the generalized defect, we argued that even in the 
presence of 1 x 1 Kronecker factors the total bound (167|) . obtained by substituting f l58| 
159]) into (|57|) . will be equal to the value ( )67|) obtained for the Kronecker product with 
absorbed lxl factors (this can be also guessed purely by analyzing ( 1671) as a function of 
sequence of sizes (rii, ..,n r )). Let us leave it as a comment that ( 1671) must, in the presence 
of 1 x 1 factors, lead to (|68|) or (1691 ) deprived of its factors (2n; — 1) corresponding to 
n\ — 1, but since then 1n\ — 1 = 1, in Theorem 13.41 one does not need to assume that 
all rife > 1. 



Theorem 3.4 The lower bound (67\ ) on the generalized defect (definition lf6o^) ) of U 
being the Kronecker product UM) with factors of size n± x n\, n r x n r 

respectively, where all > 1, is equal to 
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a) 

II ( 2 -n« - !)• ( 68 ) 

Ze{l..r} 

z/ a// rifc > 2, 

b) 

J] (2-n,-l)J . 2« 2 < 1 " r »- 1 (2^ 1 " r » + 1) , (69) 

^€{l..r}, n. ; >2 / 

if there is at least one 2 m t/ie sequence (rii, ...,n r ), where ^{L.r} is i/ie number of 
2 7 s m i/ws sequence. 

Proof 

We need to show that there holds an equality between two values, one being a relatively 
simple polynomial expression: 



a) 



n *« • n ^ - !) + !) n ^ + ^ - *)) ( 7 °) 

Je{l..r} / \l6{l..r} / Ze{l..r} 



if all nfc > 2, 

b) 

n n ] [ n ((^-!) + !) ) - ( n ^ +^-i)))-2^»- i (2^» + 1) 

v iG{l..r} / \«6{l..r} / \ie{l..r}, n,>2 J 

(71) 

if Tijfc = 2 for some k G {l,...,r}. 

and the other being the long expression in the most outer bracket in ( 16?]) . This last 
expression will be simply referred to as ( 167]) '. We will compare the above mentioned two 
quantities component by component. 

We start from the p-th component of (1671) '. by which we mean: 

e n * - 2^ i ->\^»- 1 ) • n (n,-i) 

fci<...<fc p €{l..r} \ ye{l..r}\{fci..fc p } / /e{l..r}\{fci..fc p } 

(72) 

Step 1 The, call it, left part being the result of taking only the left product from the 
inner bracket in ( 1721) is equal to: 



e n ^ n ("i - !) ) ■ ( 73 ) 

fci<...<fc p e{l..r} \/G{l..r} l€{l..r}\{ki..k p } 
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The above sum can also be found in the expansion of the left product of (I7D|) or (1711) . 
Any component of the above sum is the result of choosing, in the right subproduct 
of the left product of (170] ITT|) . p brackets at positions k\, k p from which l's will 
be taken to be multiplied with (ni — l)'s from the remaining r — p brackets. Thus 
every left part ( 1731) for p — l..(r — 1) has its counterpart in the left product of (1701 
ED). 

On the other hand, in the expansion of this product there are two components, 
namely {n\ ■ . . . ■ n r )(ni — 1) • . . . • (n r — 1) and {n\ ■ . . . ■ n r ) so far unexplained. We 
easily see that the first one is the left part of the 0-th component of (16"7]T : 

H ni - 2«^»- 1 ] • J] (n,-l) (74) 
K le{i..r} J ie{i..r} 

The second one still have to be found. 

Step 2 The, call it, right part of (172]) looks like this: 

- £ L 1 ...n fcp -2(^ 1 '-^--M)-i. "Q (ni _ 1} j (75) 

fc 1 <...<fc p g{l..r} y ie{l..r-}\{fci..fc p } / 

The expression under the sum, given for some ki, k p , can be further transformed 
into: 



II n J • 2<M 1 -»- 1 . J] ("«-!) ( 76 ) 

Jg: n fcl3 ^2} / ie{l..r}\{fei..fc p }. «i^2 

if there are no 2's in (m, ...,n r ) or if there are 2's in (n^, ...,ni r _ p ) where 
{/i,..,/ r _ p } = {l,..,r} \ {/ex,.. ,fcp}. 



II n*. U^)- 1 . J] (m-l) '2 (77) 

Jg: n fc ^2} / \ie{l..r}\{fc 1 ..fc p }, n^2 / 

if there are no 2's in (n^, ...,ni r _ p ) and there are some in (n^, ..,71^), where 
{h, ■■, l r -p} is as above. Note that in this situation the additional 2 at the end 
is caused by the fact that ^{^i, ••, U-p} = 1 due to the definition of fo. 

So, in the sum ( 1751) we add expressions like ( J76l) or ( 1771) depending on how 2's are 
scattered among the considered subsequences of (ni, ...,n r ). And we do this over 
all p = l...(r — 1) in the total sum in ( |67j) '. 

Next we will add only those summands in the sum of all right parts (175]) over 
p = l...(r — 1), which are associated with sequences (k\, k p ) containing a given 
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fixed subsequence (k qi , k qt ) = (si,...,s t ) such that n Sl , n St are all greater 
than 2, and which (these (k±, .., fc p )'s) are such that for any kj not belonging to this 
subsequence n^. = 2. The sizes n Sl , n St play the role rik q of ( !76|) and ( !77|) . All 
sequences are assumed to be increasing ones. 

• For p = t we add nothing to build (ki,...,k p ) out of (k qi , k qt ), therefore 
we have only one such summand, here q\ = l,...,qt = p. The number of 
summands is equal to (^g"'^)- Formula (176|) is used for a summand. 

• For p = t + 1, the number of summands is equal to if there are no 2's in 
(ni, n r ), or (" 2 ^ " r ^), that is the number of ways in which we can extend 
(k qi , ■■■,k qt ) to (k\, ...,k p ) by such k that n& = 2. For a summand the formula 
([76]) is used if there are at least two 2's in (m, .., n r ), that is if there is at least 
one 2 outside (n^, ...,nk p ). Otherwise formula (ITT)) is used for a summand. 

• For p = t + 2 the number of summands is if there is only one 2 in (n\, ...,n r ), 



or ( 2 ) if there are at least two. If there are more than two 2's in (ni, 
formula (jTHJ) is used for a summand, otherwise we use flTTI) 



n r 



• If p — t + ^{l-.r} and p < r — 1 the number of summands is either or 
(H}"1). In this case there are no more 2's in (ni, ...,n r ) outside (n^, ...,nk p ), 
therefore we use formula f !77p for this single summand. (This all applies if 
there are 2's in (rii, ...,n r ). If there are no 2's, we end with p = t.) 

• If p — r — 1 and ^{l-.r} = r — t, then p — t + ^{l-.r} — 1, the number of 
summands is either zero or it is equal to Vj-i) ■ Since there is one 2 left in 
(ni, ...,n r ) outside (n^, ...,nk p ) we use f |T6|) for all summands considered here. 
(This again applies if there are 2's in (ni, ...,n r ).) 

Having calculated the number of summands in each particular case, the sum of all 
of them, with sign — inherited from f |T5|) : 

-E E (n fcl ...n fep .2(fcM\^"M)-i. "Q {ni _ K 

P=l (Jfei,...,ifep) S.t. V i6{l..r}\{fci..fcp} 

.., fcqj = (si, St) 
for some q±, ..,q t 

takes the following forms. 

• If there are no 2's in the sequence (ni, ..,n r ) 

nO^^)- 1 -] n («i - 1) i ■ i . (79) 

^=1 / \«e{l..r}\{si..s t }, n ; ^2 



{7i 
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(I75|) equals to 

ik) n ^-^ • ( 8 °) 

J=l J le{l..r}\{ Sl ..s t } 

In fact, in this case it is better to analyze directly ( 175]) . where, by taking 
p — 1, ...,r — 1, one finds that any summand has its counterpart in the right 
product of (1701) . Namely, (1801) for (si,...,s t ) = (ki,...,k p ) is the result of 
choosing n\ from brackets at positions ki, k p in the right product of (1701) . 
to be multiplied by (rfy — 1) from the remaining brackets. The additional 
summands in an expansion of the right product of (T70|) are the product of all 
rifc's, ni • . . . ■ n r , which compensates for an identical product in an expansion 
of the left product of (1701) . we searched for this earlier at Step 1, and finally 
the product of all (n& — l)'s which in turn can be found in the O'th component 
(174)) of (167]) '. on the other side of the equality between (1701) and ( )6"7j) ' being 
proved. Thus we are completely done with the case when there are no 2's in 
(ni, ...,n r ). 

If the are some 2's in the sequence (ni, n r ): 

Case 1 If t + hi 1 , -, r } < r - 1 then (JTJ) equals to 

. 2 (tt 2 {l-r})-l . "Q (nj _ ^ 

ie{l..r}\{si..s t }, 

wi o r} ) - rr r> ) + ( wi 2 r> ) + + (as 

In this case we can extend (si,...,s t ) to (k±, k p ) in such a way that 
(rifej, contains all 2's from (m, ra r ). The summand associated 
with such (fci, k p ) is expressed with the use of (J77J), in which we have 
an additional factor 2. That is why expression of type (1761) standing in 
the top part of (15T1) has to be added twice, and hence the last 1 in the 
rightmost bracket of (IHTj) . This expression finally takes the form: 

-(f[n s )-2^- r ^.( J] (n i -l))-(2(Mi-}) + i) . 

\i=l / y6{l..r}\{si..s t }, "i^2 / 

(82) 

Case 2 If t + hi 1 , ■-, r ) = r , that is t = r - fl 2 {l, r} we get for (1781 

a=l / 2G{l..r}\{si..s t }, r H +2 
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In other words, outside the subsequence (n Sl , n St ) there are only 2's 
in the sequence {n\, n r ). Every summand associated with (ki,...,k p ) 
produced out of such (si,...,s t ) is expressed with the use of (!76|) . Note 
that in the case of only one 2 in (m, n r ) the sum of binomials containes 
only (q) . The sum ( 183]) takes the final form: 

(fe{i..r})-i . I Y[ (n,-l) J • (2( fe ^-^-l) = 

«6{l..r}\{si..s t }, / 

) 2^»- 1 • i n (m - 1) ] ■ (2 (wi - r}) + 1) + 

v*=l / y6{l..r}\{si..s t }, n ( ^2 / 

]]^)-2^). ( H (n,-l))(84) 

v*=l / \(e{l..r}\{ Sl ..s t }, n ; ^2 / 

The reader should realize, that the Case 2 occures only once, when (si, St) 
is of maximal length. Then we get an additional component, the second part 
of ( 184|) . In general we add expressions like ( 182|) . over all possible nonempty 
sequences (si, s t ), that is over all t ranging from 1 to r — ($ 2 {l..r}). So let 
us first concentrate on them. 

As in the case without 2's in (ni,...,n r ) we note that any expression like 
(182]) has its counterpart in the right product of (ITT]) . Again it is a matter of 
choosing brackets from which we take n^'s, at the positions corresponding to 
the positions of si, s t in the increasing subsequence (k : > 2) of (1, r). 
Note that in this process we also deal with the situation when all > 2 are 
used in (182]) . in fact in the 'standard' part of (184]) ((s l7 s 4 ) of maximum 
length, Case 2 above), which corresponds to taking only n^s from brackets 
in (J7TJ. 

What remains in an expansion of ( ITT]) are the values: 



(n, - 1) • 2« a * 1 " r »- 1 (2«^ 1 " r » + 1) and (85) 

^G{l..r}, n ; >2 y 

n n <> ( 8e ) 

Z€{l..r} 

where the last product was mentioned as unexplained at the end of Step 1. 
On the other side, in (l6Tj) '. we have still left the right part of the 0-th compo- 
nent (IT4]) of (I67p '. which can be written as: 

J] (ni-l)\ ■2^ 1 - r ^- 1 , (87) 

ie{l..r}, ri[>2 I 
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the additional component in (I84p occuring only for (si,...,s t ) of maximum 
length, therefore written as 



J] m 2 {h{1 ~ r}) , (88) 

ni>2} j 

and finally all summands in the left parts ( 175|) for p = 1, r — 1 which are as- 
sociated with sequences (ki, k p ) such that there are only 2's in (n kl , ...,n kp ). 
We add these summands in Step 3 which completes the proof. 

Step 3 The expression under the sum in ( 1751) . in case of (k\, k p ) such that (n^, ...,nk p ) 
(2, 2), takes the form 

2 (fc{1 " r}) - 1 • {m - 1) (89) 

{l: n,>2} 

if the number of 2's in (jii, n r ) is greater than p, 



2 C» a {i--r»-i . JJ („, _ 1) . 2 (90) 

\{«: n ; >2} / 

if the number of 2's in (n±, ...,n r ) is equal to p, p = (^{l-.r}). This is again 
caused by the definition of jfe. In this case 2™ 2 * 1 " r }\{ kl " kp » is equal to 2 while 
there are no 2's among n k indexed by k G {l. r} \ {ki..k p }. 

Now we calculate how many summands of this type can be found in the sum of all 
the right parts (1751) . It is assumed that at least one 2 can be found in (m, ■■■,n r ). 
The other case has been completely resolved in the paragraph following formula 

USD)- 

• For p = 1 the number of summands is equal to (" 2 ^" r ^), the number of ways 
we choose a one element subsequence, containing only 2's, from (ni,...,n r ). 
We use (19"0|) for a summand if there is only one 2 in {n X) ■■■,n r ), otherwise we 
use 



For p = 2 the number of summands is either if there is only one 2 in 
(ni, n r ), or (^I"^)) we choose a two element subsequence this time. We 
use (190]) for each summand if there are only two 2's in (m, ...,n r ), or we use 
(1591 if there are at least three. 



For p = $2{l..r} the number of summands is (^{i r}) anc ^ we °f course use 



formula (190|) for this single summand. 
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We add all the considered summands below, where the negative sign is inherited 
from (!75|) . In the first sum we still use the form of a summand originating from 
73): 



r-l 



E n fcl ...^.2^-^-M)-i. Y[ (m-1) 

P=1 (k 1: k p ) S.t. V «6{l..r}\{ fcl ..M 

(rifc^ nfc ) = (2, 2, 2) 

-•■--'■n >i (-.){( Hi 1 jl ) + ( wi 2 - r> ) + - - ({:{";!) + (9i) 

where the last 1 is the result of applying formula ( 1901) . with that additional mul- 
tiplication by 2, to the single summand for p = ^{L.r}. The final result is 
thus: 

J] {m - 1) ] • 2«^ 1 - r »- 1 • 2« 2 ^- r » (92) 

n ; >2} / 

This is the end. What is left for us to say is that: 



92]) together with (1H7|) . both within ( )67|) '. compensate for (|85|) being part of 



88|) . a part of (1671)'. is equal to fl86|) coming from ((71 



Our calculation of the bound fl^S j loT?]) on the generalized defect D(f/) of t/ given by ffT2"|) 
is based on the formula for the dimension of Mu given by Theorem 13.21 which in turn 
is based on the decomposition of M.jj into direct sum components provided by Theorem 
12.61 As far as the case with no 2's in the sequence of sizes (rai, ...,n r ) is concerned, we 
can use a simpler reasoning leading to the corresponding bound (|68|) . To this end let us 
consider the following lemma, in which we use spaces defined in ( |66l) . 



Lemma 3.5 Let U and V be unitary matrices of size N x N and M x M respectively. 
Let R G \fu and S G Vy, where R, S are real matrices of the size identical with that of 
U, V, respectively. 

ThenR®S G Vu®v- 



Proof 

From the assumptions there holds 

iRoU = EU and iS o V = FV (93) 
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for some antihermitian matrices E, F. Then 

(iRoU)®(iSoV) = -(R®S)o(U®V) = EU®FV = (E®F)(U®V) , (94) 
from which we get 

i(R®S)o(U®V) = (-iE®F)(U®V) , (95) 
where (— iE ® F) is antihermitian because E <g> F is hermitian. So R <8> 5* G Vy^y. ■ 
Therefore, if \R) are bases for Vn-r*o, = l..r, then R) ® . . .® R) form a 

V 4 A=1..D([/W) u ' »i 

set of D(U^) ■ . . . ■ D(U^) independent vectors within Vf/a)®...®^). Hence: 

Corollary 3.6 For a Kronecker product of unitaries ( fl^jj the generalized defect is super- 
multiplicative: 

D (If® <g> . . . ® C/ (r) ) > D([/«) • ... -D(C/W) (96) 

The above mentioned set of i?^ (g> . . . <S> -R^ needs not to be a basis for Vyci^.^cr). 
For example, for the 2x2 unitary Fourier matrix: 

we have that D(F 2 <g> F 2 ) > D(F 2 ) • D(F 2 ), as D(F 2 <g> F 2 ) = 10 and D(F 2 ) =3. 

The space V^so constructed for a unitary C/( fc ) of size rik x contains at least 2n^ — 1 
independent real matrices, denoted further by r\ , i = l..(2n& — 1), which are all zero 
matrices except for the /-th row or m-th column filled all with ones, where I = l..rif. and 
m = 2..rik- The corresponding directions iR^ o were mentioned at ([2]) . Thus 

D(t/«) > 2n k -l , (98) 

and consequently, by Corollary 13. 6[ 

D (U {1) (8) ... (8) U {r) ) > (2m - 1) ■ (2n r - 1) . (99) 

This is equal to (16"8|) . which is also valid if there is only one 2 in the sequence of sizes 
(ni, ...,n r ). 

If there is more than one matrix of size 2x2 among U<®, in (|99D we get a lower bound 
which is smaller than the lower bound ( 1691) . This is caused by the fact that, for x denoting 
the number of 2's in (rii, ...,n r ), there holds an equality for natural x > 1: 

2 x-i . + y > (2 . 2 - l) x , (100) 

which is equivalent to 

4 X -3 X > T-2 X , (101) 



1 1 
1 -1 



(97) 
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which in turn is true because, after dividing both sides by (4 — 3) = (3 — 2) = 1, on 
the left the summands are greater then on the right. So, in this situation, unsplitting the 
2x2 matrices one by one from . . . <g) is not a good strategy and it is better 
to write (by Corollary 13.61 and by ( J98l) ): 



D (TJM ® . . . ® Z/M) > ( J] B(U^)) D ( (g)f/( fc M > 

\n k >2 J \n k =2 J 

Y[(2n k -l)).D((g)uA . 

n k >2 J \n k =2 J 



(102) 



It would be interesting to conceive a set of 2 X 1 (2 X + 1) independent matrices of some 
elegant structure within _{/(*), where x is the number of 2's in (ni, ...,n r ), so that 
we could write (we can, but basing on Theorem 13.41) : 




> 



IX— 1 



(2 X 



(103) 



(We mean here a guess similar to the above choice of (2rik — 1) matrices in V^j 
This, together with H102j) . would provide the bound ( 1691 ) of Theorem 13.41 probably more 
directly than with the use of Theorems 13.21 and 12.61 Definitely, taking the 3 X matrices 

<8> 



R?f* , where factors come from the set 

n k =2 i k ' 



1 1 







1 1 



1 
1 



(104) 



is, due to (HOOp . not enough. 



4 Conclusions 

In this article we have achieved two things. 

The first seems to be of great importance in numerical calculation of the defect (Definition 
II. ip for large unitary NxN matrices U with a Kronecker product structure. This involves 
calculation of the dimension of a certain space Mjy (defined in Q), effectively in H N , 
spanned by the set A4.u (defined in (llOp). associated with U. It amounts to calculating 
the column rank of a certain matrix, call it M(U), built on the entries of U, and can 
be performed for example by using rank/svd of MATLAB. When U has a Kronecker 
product structure, M;y can be split into a direct sum with a large number of components 
(see Theorem 12. 6p . So, instead of applying the procedure to M(U) we can apply it to a 
number of its submatrices. This will make the whole calculation more reliable and may 
prevent divergence in svd used in calculation of the rank. 
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The second thing we have achieved in this work is the calculation, based on the above 
result, of the lower bound on the generalized defect D(f) (defined in (1651) . being the defect 
plus 2N—1) when U is a Kronecker product. The generalized defect is supermultiplicative 
with respect to Kronecker subproducts of U (Corollary 13. 6 1) and this allows us to trivially 
retrieve the lower bound of Theorem 13 .41 when the number of 2 x 2 Kronecker factors does 
not exceed 1. In the other case Theorem 13.41 gives a better bound than that obtained 
using supermultiplicativity of D(U). We conjecture that, in either case, the lower bound 
is attained by most matrices with a fixed Kronecker product structure (fixed, up to an 
order, sequence of sizes of Kronecker factors). 

All the formulas, the one expressing the direct sum forming My (see (]52p ). the resulting 
one expressing the dimension of My (see (1571) ). and the one expressing the lower bound 
(see (1671) . in a compact form (168|) or (1691) ) are valid also when there are lxl Kronecker 
factors among those forming U, yielding the correct values we would get if the lxl factors 
were absorbed into larger factors. This property allows, in calculation of dim (My) or its 
bound, the use of a procedure requiring a fixed number of factors, where the lxl [l]'s 
can be taken as potentially missing factors to extend a shorter Kronecker product. 

The author believes that the first of the above mentioned achievements will open the 
way to easier calculation of the defect for large matrices U with a Kronecker product 
structure. Such calculations may be performed in order to asess whether a given U gives 
rise to a smooth family of inequivalent unitaries (i.e. not obtained one from another by 
multiplying rows and columns by unimodular numbers) with the moduli of entries fixed 
at the values sitting in [iC/ij]^ - =1 N , and how large its dimension can be. This is a part 

of a question about unitary preimagies V to a doubly stochastic matrix [|C/j,j| ] 4 ._ a N , 
that is about V such that \Vi j\ = |£/y| for all This motivation is wider described 
in a prequel [I] to this paper. A special case, [I^M'^Lj-i n = [V-^1ij=i jy ^ ne searcn 
for complex Hadamard matrices, attracts even more interest, and it was dicussed in the 
Introduction. 
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